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1.  INTRODUCTION 


There  are  three  basic  versions  of  the  finite  element  method,  called  the  h,  p 
and  h-p  versions.  They  are  essentially  characterized  by  the  way  in  which  the  finite 
element  meshes  and  polynomial  degree  of  elements  are  chosen.  They  differ  in  computer 
Implementation  (program  architecture)  and  mathematical  analysis.  This  paper  is 
concerned  mainly  with  the  question  of  how  the  meshes  and  polynomial  degree  of  the 
elements  affect  the  accuracy  of  finite  element  solutions.  Our  approach  is  to  fix 
certain  parameters  or  their  relation  and  increase  the  number  of  degrees  of  freedom  so 
that  the  finite  element  solutions  converge  to  the  exact  solution.  Such  a  systematic 
increase  of  the  number  of  degrees  of  freedom  is  called  extension  because  it  can  be 
interpreted  as  a  systematic  extension  of  finite  element  spaces.  g.,- 

When  emphasis  is  on  analysis  of  accuracy  and  not  aspects  of  implementation,  then 
we  speak  about  the  h,  p  and  h-p  extensions  rather  than  versions.  Understanding 
the  various  extension  processes  and  their  numerical  performance  is  essential  for 
resolving  certain  basic  questions  of  implementation. 

The  h-extension  is  the  most  commonly  used  approach  to  error  reduction.  The 
polynomial  degree  (p)  of  the  elements  is  fixed  and  the  errors  of  approximation  are 
reduced  through  mesh  refinement.  The  size  of  the  elements  is  usually  denoted  by 
h,  hence  the  name:  h-extension.  Typically,  the  polynomial  degree  of  elements  is 
low,  usually  p  -  1  or  p  =*  2. 

In  the  p-extenslon  the  mesh  is  fixed  and  convergence  is  achieved  by  increasing 
the  polynomial  degree  of  elements  either  uniformly  or  selectively. 

The  h-p  extension  combines  the-  h-  and  p-versions,  i.e.,  reduction  of  error  is 
achieved  by  mesh  refinement  and  concurrent  choices  in  the  polynomial  degree  of 
elements. 

The  parameters  that  characterize  extension  processes  can  be  chosen  either _a 
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priori ,  on  the  basis  of  certain  characteristics  of  the  exact  solution,  known  a 
priori,  or  a  posteriori,  through  utilization  of  some  feedback  procedure  in  which  case 
the  parameters  of  the  extension  process  depend  on  previously  computed  data. 

The  analysis  and  especially  optimization  of  extension  processes  (selection  of 
optimal  meshes,  polynomial  degree  distributions,  etc.)  presented  herein  indicate  the 
potential  of  alternative  approaches  and  provide  a  basis  for  decisions  concerning 
Implementation. 

In  order  to  keep  the  essential  points  in  focus,  we  consider  only  two  simple 
model  problems  based  on  the  displacement  formulation  and  two  measures  of  error:  The 
error  measured  in  energy  norm  and  the  error  of  stress  components  computed  at  specific 
points. 

Specifically,  we  denote  the  exact  and  finite  element  solutions  respectively  by 
Uq  and  u.  The  error  is  then  e  *  ug  -  u.  The  energy  norm  of  e  is  denoted  by 
lelg  and  is  defined  as  the  square  root  of  the  energy  of  thd  error: 

lel£  =  We))1^  . 


The  relative  error  in  energy  norm  is  denoted  by  lei  and  is  defined  as: 


lel 


ER 


lelE 

'"o'e 


The  error  in  stress  components  at  some  poin  x0  is  defined  as 


eij<xo>  =  l^v  -  vv1 

where  °ij(xo)  and  respectively  denote  the  exact  and  computed  components 

of  the  stress  tensor  at  point  Xq.  The  relative  error  in  stresses  is  defined  by: 


e 


R 

ij 
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leij  (XQ>I 
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The  one  dimensional  problem  can  be  analyzed  theoretically  and  experimentally  in  great 
detail.  One  dimensional  problems  can  also  serve  as  models  for  higher  dimensional 
problems  which  are  vastly  more  complicated  and  less  well  understood.  Presentation  of 
details  and  derivation  of  formulas  quoted  herein  is  beyond  the  scope  of  this  paper. 
For  the  proofs  of  theorems  in  one  dimension,  we  refer  to  [1],  [2j.  For  further 
details  and  application  to  two  dimensions,  we  refer  to  [3],  [4],  [ 5 ] ,  [6],  [8J. 


2.  MODEL  PROBLEMS 


2.1.  Model  Problems 


We  consider  the  following  simple  model  problems: 


-u"(x)  -  f(x),  x  €  I  -  (0,1) 


(2.1) 


u(0)  -  u(l)  -  0 


(2.2) 


and  exact  solutions  of  the  form: 


uQ(x)  •  (x+£)°  -  K°  -  x[(i+oa  -  S°]  ( 

with  a  >  V2  and  £  >  0.  The  solutions  minimize  the  potential  energy  defined  as: 


(2.3) 


x(u)  =  W(u)  -2  /  fu  dx 

0 

1 

where  W(u)  *  J  (u')2  dx. 

0 

The  finite  element  solutions  are  characterized  by  the  mesh  and  p- 
distributions.  Specifically  we  denote  the  mesh  by  the  partition: 


A  =:  0  = 


A  .  A  »  y  A 
x0  <  xl  <  —  <  XM(A) 


where  xA,xA  ’***  ,XM(A)  are  t*ie  mes^  (°r  n°dal)  points.  The  j1"*1  finite  element  is 

denoted  by:  Ij+j  =  (xj»xj+i)*  The  size  of  the  j1'*'  element  is  defined  by  hj 

»  xj  -  The  size  of  the  largest  element  is  denoted  h(A).  The  set  of  all 

functions  w,  defined  on  I  that  satisfy  the  following  conditions: 

(a)  W(w)  <  •  , 

(b)  the  boundary  conditions  (2.2), 

(c)  on  1^  w  is  a  polynomial  of  degree  pA, 

P  ( A)  a  a  A 

is  denoted  by  S  (A)  where  p(A)  -  fp.  ,p0, . . .  ,p„,  .  N  J  is  the  vector  of  p- 

2(A)  1  2  M(A) 

distribution.  S  (A)  is  called  the  finite  element  space.  The  number  of  degrees 

of  freedom  is  denoted  by  N(A):  r 


P(A) 

N(A)  =  dim  S  (A) 


M(A)  . 

I  p,1  - 1. 


2(A) 


The  finite  element  method  consists  of  finding  u  £  S  which  minimizes  x(u)  over 

_e(a) 

S  (A).  The  extension  processes  are  characterized  by  the  selection  of  sequences 
of  A  and  £(.&)• 
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2.2.  The  Two-Dimensional  Model  Problem 


S- 


h: 


We  shall  consider  the  model  problem  of  Fig.  2.1  assuming  plane  strain  conditions 
and  using  the  elastic  T*arameter9  E  ■  1,  v  ■  0.3.  The  tractions  on  the  boundary  are 
chosen  so  that  the  exact  solution  is  known  and  the  singularity  at  the  crack  tip  is 
characterized  by  the  stress  intensity  factors  KI  “  KII  “  l* 

The  solution  of  Model  Problem  1  has  certain  similarities  with  the  solution  of 
Model  Problem  2  from  the  point  of  view  of  numerical  performance  when  {  *  0  and  0.5 
<  a  <  1. 


3.  THE  PERFORMANCE  OF  THE  FINITE  ELEMENT  METHOD  IN  ONE 
DIMENSION  WITH  RESPECT  TO  THE  ENERGY  NORM  MEASURE 


We  now  demonstrate  how  the  performance  of  the  various  extension  processes 
depends  on  the  mesh  and  the  polynomial  degree  of  elements.  The  principal  mathematical 
tool  is  asymptotic  analysis  which  provides  information  on  how  the  error  depends  on 
the  mumber  of  degress  of  freedom  when  only  one  parameter  is  being  varied,  and  on 
certain  essential  characteristics  of  the  exact  solution  provided  that  the  number  of 
degrees  of  freedom  is  large.  We  shall  consider  Model  Problem  1  and  choose  the 
parameters  £  and  a  so  as  to  represent  problems  with  solutions  of  various 
smoothness  (£  ■  0,  a  >  V2  )  and  £  >  0. 


3.1.  H-Extensions  Based  on  Uniform  Mesh  Refinement  (£  *  0) 

In  this  case  M(A)  *  p  is  fixed.  The  number  of  elements  fully  charac¬ 

terizes  the  mesh  and  the  error  in  energy  norm  for  £  “  0  is  estimated  as  follows: 


|ei  t 

ER  MV 

(3.1) 

where  the  symbol  *  means  "asymptotically  equal",  C(a) 
of  the  mesh  and  the  polynomial  degree  of  elements: 

is  a  constant  Independent 

p  -  min(or-  V2  »p) 

(3.2) 

p  ■  2a  -  1 

(3.3) 

and. the  number  of  degrees  of  freedom  is  given  by: 

N  -  Mp-1 . 

(3.4) 

Although  (3.1)  is  an  asymptotic  estimate,  it  holds  even  for  reasonably  low  values 
of  N.  In  order  to  demonstrate  this,  we  have  computed  the  value 


D  =  lei ER  Mwpp.  (3.5) 

D  is  called  the  numerical  value  of  C(a).  The  results  for  a  »  0.7  (and  therefore 
p  -  0.2,  p  *  0.4)  are  shown  in  Table  3.1.  *e*ER  *s  s^own  as  "percent  relative 
error".  It  is  seen  that  (3.1)  holds  well  also  when  M(A)  is  small  and  the  error  is 


Li' .  i 


Vo 


large.  In  view  of  the  fact  that  p  Is  twice  the  value  of  p,  for  the  same  number 
of  degrees  of  freedom  the  higher  order  elements  perform  better. 

3.2.  H-Extensions  Based  on  Nearly  Optimal  Mesh  Refinement  (£  *  0) 

We  once  again  consider  the  case  £  ■  0.  In  this  case  the  optimal  mesh  (for 
fixed  p)  is  asymptotically: 


where: 


ffc*B 


p  +  *4 


(3.6) 


(3.7) 


The  estimate  of  error  measured  in  energy  norm  is  given  by: 


P 


'1_M-[(2a-l)g-2p] 

(2a-l)B-2p 


=  C  N'P 


provided  that 


B  >  B  = 


(3.8a) 


(3.8b) 


When  B  <  B  then  the  rate  of  convergence  with  respect  to  M  decreases.  The 
reliability  of  estimate  (3.8)  for  B  =  Bq  is  shown  in  Tables  3.2a  and  3.2b. 

With  reference  to  Tables  3.2a,  3.2b,  we  note  the  following  observations: 

(a)  The  asymptotic  estimate  (3.8)  is  of  good  quality  when  M  >  2p.  For  M  < 

2p  the  formula  is  pessimistic.  The  reason  for  this  is  that  (3.8)  is  based  on  the 
assumption  that  M  -*•  «>.  Therefore  it  cannot  be  expected  to  give  close  estimates  for 
low  values  of  M.  In  the  case  M  <  2p,  analysis  based  on  p-extension  rather  than 
h-extension  should  be  used. 

Table  3.3  shows  the  error  for  optimal  distribution  of  the  nodal  points  for  M  ■ 

2  and  a  =  1.1.  It  is  seen  that  for  small  values  of  p  (p  *  1,2)  the 
asymptotically  optimal  mesh  performs  very  nearly  as  well  as  the  optimal  mesh. 

(b)  When  a  is  small  (strong  stress  singularity  occurs)  the  optimal  refinement 
is  so  strong  that  roundoff  limitations  are  encountered  even  when  the_computations  are 
performed  in  double  precision.  Table  3.4  shows  the  values  of  Bq»  B,  the 

coordinate  of  the  first  nodal  point  for  the  optimal  mesh,  x^  (opt),  and  for  the  mesh 

at  which  the  rate  of  convergence  begins  to  decrease,  x^  (min),  for  a  *  0.7,  M  =>32. 


(c)  Overrefinement  is  more  advantageous  than  underrefinement.  If  a  is  not 
known  precisely,  then  the  refinement  should  be  designed  for  lower  bound  estimates  of 
a.  Over refinement  increases  the  value  of  2  in  (3.8a)  but  does  not  alter  the  rate 
of  convergence  (N~P).  The  penalty,  in  terms  of  increased  values  of  C  and 
increased  values  of  N  for  achieving  comparable  levels  of  accuracy  are  shown  in 


SCOT 


*■»  •«»•» 

K* .•  «*  ■  ' 


Table  3.5.  The  mesh  refinement  is  optimal  for  a  -  0.7,  the  penalty  values  are 
shown  for  a  *1.1,  1.6,  2.1. 


It  is  seen  that  the  penalty  for  over refinement  is  not  small.  Underrefinement 
also  reduces  the  rate  of  convergence.  Therefore  the  penalty  for  underrefinement  is 
still  larger.  This  shows  the  importance  of  correct  selection  of  mesh  refinement. 


3.3.  P-Extensions  Based  on  Uniform  Meshes  (£  “  0) 

In  this  case  M(A)  «  p  and  the  estimate  (3.1)  holds.  The  results  of 
computational  experiments  are  shown  in  Table  3.6. 

On  comparing  the  results  on  the  basis  of  the  number  of  degrees  of  freedom  (N  * 
Mp-1),  it  is  seen  that  the  best  choice  is  M  ■  1.  Comparing  Tables  3.6  and  3.2a  it 
is  evident  that  h-extensions  based  on  optimal  meshes  yield  better  results  than  p- 
extensions  based  on  uniform  meshes.  The  performance  of  p-extensions  cannot  be 
improved  substantially  through  optimizing  the  p-distribution.  Table  3.7  shows  the 
effect  of  optimal  p-distributions  for  M  *  2  and  a  *  0.7,  1.1. 


3.4.  H-p  Extensions  (£  *  0) 

In  this  case  the  meshes  and  p  distributions  are  optimized  concurrently.  The 
assymptotically  optimal  mesh  is  characterized  by  the  following  geometric  progression 
of  nodal  points: 

xj  ■  ko<A)  J  j  “  (3.9) 

where: 

kq  -  (/I  -  l)2.  (3.10) 

The  polynomial  degrees  of  elements  are  assigned  as  follows: 


Pj  *  [ s(a) j ]  (3.11) 

where 

s(a)  -  2(a-V2)  (3.12) 

and  [•]  means:  the  integral  part. 

The  error  estimate  is: 

lelER  -  C(a)[(/2  +  l)2]-^°-  1/2  )N 

-YQ/(a-  V2)N 

-  C(a)e  ,  Y0  - 


1.574 


(3.13) 
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The  numerical  performance  of  h-p  extensions  is  illustrated  in  Table  3.8.  On 
comparing  the  results  with  those  in  Table  3.2a(  it  is  seen  that  the  rate  of 
convergence  of  the  h-p  extension  is  much  greater  than  the  rate  of  convergence  of 
the  h-extension  based  on  asymptotically  optimal  meshes. 


3.5.  The’  h-p  Extension  with  Uniform  p  and  Optimal  Mesh  (£  -  0) 

The  mesh  is  the  same  as  in  Section  3.4,  i.e.,  with  the  nodal  points  defined  by 
(3.9),  (3.10)  and  we  assume  that  the  p-distribution  is  uniform.  In  this  case  the 
estimate  for  <  a  <  1  is: 


. . . 

/ (a-  V2  )N 

,e,ER  -  C(a)-2^i'/2  (3-U> 

P 

where,  as  in  (3.13),  ■  1.574  and: 

p  -  [  s(a) ]M(A)  (3.15) 

s(a)  -  2(o  -V2).  (3.16) 


The  numerical  performance  of  the  h-p  extension  with  uniform  p  is  shown  in  Table 
3.9.  On  comparing  Table  3.9  with  Table  3.8,  it  is  seen  that  the  performance  of  the 
h-p  extension  with  uniform  p  is  not  substantially  different  from  that  of  the  h-p 
extension  with  optimal  p-distribution.  The  performance  can  be  analyzed  also  for  p- 
distributions  other  than  that  given  in  (3.15).  When  p  increases  more  rapidly  than 
(3.15),  then  the  rate  of  convergence  diminishes  until  it  reaches  the  algebraic  rate 
characteristic  of  p-extensions.  When  p  increases  less  than  (3.15),  then  the  rate 
of  convergence  diminishes  because  h(A)  does  not  change. 


3.6.  H-Extensions  that  Utilize  Feedback  (£  ■  0) 

In  Section  3.2  it  was  pointed  out  that  the  quality  of  performance  of  h- 
extensions  depends  on  the  mesh  design.  Proper  mesh  design  depends  on  the  exact 
solution  which  generally  is  not  kngwn.  It  is  possible  to  devise  feedback  procedures 
however,  that  construct  meshes  which  asymptotically  perform  as  well  as  the  optimal 
meshes.  Such  feedback  procedures  are  called  adaptive  [9],  [10],  [11].  [12]. 

Tables  3.10a,  3.10b  show  the  results  of  numerical  experiments.  The  numerical 
value  D  is  based  on  a  formula  for  optimal  meshes  that  utilize  only  nodal  point 
which  can  be  constructed  by  successive  bisection  of  elements,  not  all  meshes  as 
considered  before,  because  the  feedback  procedure  uses  only  such  meshes. 


3.7.  H-p  Extensions  that  Utilize  Feedback 

It  is  possible  to  devise  feedback  procedures  that  perform  nearly  as  well  as  the 
optimal  h-p  extension.  Results  obtained  with  such  a  procedure  are  shown  in  Table 
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3.11  (o  -  0.7).  The  meshes  were  generated  by  bisection,  therefore,  D  is  based  on 
an  estimate  developed  for  such  meshes  only. 


3.8.  p-Extension  Based  on  Properly  Designed  Meshes  and  Feedback  (5  *  0) 

In  this  case  the  mesh  is  strongly  graded  toward  singular  points,  on  the  basis  of 

(3.9),  with  iCq  ■  (/2  -  1;  *  0.1715  or  slightly  smaller,  say  Kq  -  0.15  and  M(A) 
is  fixed.  Theupolynomial  degree  of  elements  is  uniformly  increased.  The  error 
decrease  first  exponentially  and  if  p  is  too  large,  then  algebraically,  as 
explained  in  Sec-tion  3.5.  Feedback  is  utilized  to  ensure  (through  proper  selection 
of  M(A))  that  the  desired  accuracy  is  reached  in  the  range  where  the  convergence  is 
exponential . 


3.9.  Smooth  Solutions  (£  >  0) 


We  have  considered  various  extension  processes  when  the  solution  has  singular 
character.  When  the  solution  is  smooth,  then  p-extensions  perform  especially  well 
for  small  M.  The  error  estimates  for  5  >  0  in  (2.3)  are  as  follows: 


(a)  for  C  -  0 


(b)  for  £  >  0 


where: 


lei 


ER 


C(a) 

2o-l 

P 


lei 


.3^ 


'T+c  -  SI 
T h  +  /e 


(3.17) 


(3.18) 


(3.19) 


In  Table  2.11  results  are  presented  for  a  =*  0.7  and  £  =  0,  £  *  0.01  and  £  ■ 

0.1.  These  results  demonstrate  that  the  performance  of  p-extensions  very  rapidly 
improves  with  increasing  smoothness  of  the  solution. 


A.  ACCURACY  OF  STRESS  APPROXIMATIONS  IN  ONE  DIMENSION 


In  one  dimension  the  stress  is  simply  u'(xq).  In  contrast  to  the  two 
dimensional  case,  the  behavior  of  the  finite  element  solution  in  our  example  is 
entirely  local,  therefore  we  need  to  consider  only  the  case  with  one  element.  The 
results  of  numerical  experiments  for  a  =  0.7  and  various  £  values  are  shown  in 
Table  4.1.  It  is  seen  that  the  element  that  contains  the  singularity  (5  =0)  yields 
very  poor  stress  approximations. 

When  quantities  other  than  the  energy  are  of  interest,  for  example  stresses, 
then  the  mesh  and  p-distribution  should  be  optimal  or  nearly  optimal  with  respect  to 
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85 

u 

1 


the  purpose  of  computation.  Optimal  meshes  and  p-dlstrlbutions  depend  not  only  on 
the  purpose  of  computation  but  also  on  the  method  used  for  computing  the  quantities 
of  Interest.  See,  for  example,  [10], {13]. 


5.  PERFORMANCE  OF  THE  FINITE  ELEMENT  METHOD  IN  TWO 
DIMENSIONS.  WITH  RESPECT  TO  THE  ENERGY  NORM 


The  theory  of  two  dimensional  problems  presents  more  difficulties  and  Is 
therefore  less  well  understood  than  the  theory  of  one  dimensional  problems. 
Nevertheless  there  are  Important  similarities  which  make  it  possible  to  gain  valuable 
Insight  from  the  analysis  of  one  dimensional  problems.  There  are  important 
differences  also,  for  example,  error  in  stresses  in  two  dimensions  behave  quite 
differently  from  the  errors  in  stresses  in  one  dimension.  Stress  computations  are 
discussed  in  the  next  section.  Here  we  discuss  the  properties  of  various  versions 
and  present  numerical  results  for  our  two  dimensional  model*  problem.  The  results 
were  obtained  by  means  of  the  h-version  program  FEARS  with  feedback  capabilities 
[14]  and  the  p-version  program  FIESTA- 2D  [15].  FEARS  has  elements  of  polynomial 
degree  1  only,  the  polynomial  degree  of  elements  in  FIESTA-2D  range  from  1  to  8. 


5.1.  H-Ex tens ion  Based  on  Uniform  Meshes 


The  estimate  for  our  model  problem,  defined  in  Section  2.2,  is: 


C(p)N 


(5.1) 


Detailed  theoretical  analysis,  comparable  to  the  one  dimensional  case,  is  not 
available.  The  results  of  computations  are  shown  in  Table  5.1,  where  D  represents 
the  numerical  value  of  C(p).  The  results  indicate  that  the  asymptotic  estimate 
(5.1)  is  of  good  quality. 


5.2.  H-Extension  with  Feedback 


As  in  the  one  dimensional  case,  the  sequence  of  optimally  designed  meshes  leads 
to  a  rate  of  convergence  independent  of  the  singularity.  The  estimate  for  optimally 
designed  meshes  is: 

-  C(p)N-p/2 


(5.2) 


Note  that  the  exponent  of  N  is  p/2,  not  p  as  in  the  one  dimensional  case. 


H-extension  with  properly  utilized  feedback  (adaptive  approach)  should  lead  to 
the  same  asymptotic  rate  of  convergence  as  the  optimally  designed  sequence  of 
meshes.  Table  5.2  shows  the  results  obtained  with  FEARS. 


5.3.  P-Extension  on  Uniform  Mesh 


In  this  case  the  estimate  is  [3]  [4]: 


where  e  >  0,  arbitrary.  It  is  not  known  whether  the  term  e  can  be  removed 


Table  5.3  shows  the  results,  obtained  with  FIESTA-2D  using  four  square  elements. 


5. A  H-p  Extension 

The  estimate  for  h-p  extensions  based  optimal  mesh  and  on  either  optimally  or 
uniformly  distributed  p  is: 

,e,ER  ‘  C  eV"  ' 

In  the  two  dimensional  case  the  optimal  value  of 
3 

whether  the  term  /N  can  be  improved.  The  value 
of  p. 

5.5.  p-Extension  Based  on  Properly  Designed  Mesh  and  Feedback  Information 

As  in  the  one  dimensional  case,  the  p-extension  performs  much  the  same  way  as 
the  h-p  extension  when  p  is  not  two  large  and  the  mesh  is  properly  designed.  For 
large  p  the  p-version  performs  as  if  the  mesh  were  uniform.  An  example  is 
presented  in  Table  5.4.  The  fact  that  the  rate  of  convergence  slows  for  high  p  is 
an  indication  that  the  mesh  should  be  refined.  Slowing  of  the  rate  of  convergence 
can  be  detected  which  is  the  feedback  information  needed  for  increasing  the  number  of 
finite  elements. 


T  >  0.  (5.4) 

Y  is  not  known  nor  is  it  known 
of  y  depends  on  the  distribution 


5.6.  Smooth  Solutions 

When  the  solution  is  smooth,  the  p-version  is  very  effective  and,  as  in  the  one 
dimensional  case,  the  convergence  is  exponential. 


6.  STRESS  COMPUTATIONS  IN  TWO  DIMENSIONS 


Stress  approximations  behave  quite  differently  in  two  dimensions  than  in  one 
dimension.  In  one  dimension  in  our  example  the  error  depended  only  on  the  behavior 
of  the  solution  in  the  particular  element  in  question,  i.e.,  the  error  was  completely 
localized.  In  two  dimensions,  on  the  other  hand,  the  error  is  comprised  of  two 
parts:  the  local  error  and  the  error  associated  with  all  other  elements.  This 

second  part  is  called  pollution  error. 

The  error  in  stresses  depends  to  a  large  extent  on  how  the  stresses  are 
computed.  Indirect  techniques  are  available  which  substantially  reduce  both  the 
local  and  pollution  error,  as  compared  with  the  conventional  (direct)  methods  of 
stress  computations  [7],  [10],  [13]. 


6.1.  Performance  of  h-Extensions  Based  on  Nearly  Optimal  Meshes 


Let  us  examine  the  stresses  at  point  *  x2  “  0.25  in  the  model  problem. 
Solutions  were  obtained  by  means  of  the  computer  program  FEARS.  The  point  under 
consideration  is  the  vertex  of  four  elements.  Therefore,  four  different  values  can 
be  computed,  using  the  derivatives  of  the  four  elements  and  the  appropriate  stress- 
strain  low.  The  relative  errors,  are  shown  in  Table  6.1  for  the  three  stress 
components  o^,  a22*  °12  ^or  t*ie  *our  adjacent  elements.  The  error  of  the 

average  value  (A)  is  also  shown.  It  is  seen  that  the  error  of  the  average  value  is 
smaller  than  the  error  in  most  elements.  This  is  a  well  known  fact  which  is 
generally  utilized  in  stress  computations. 

The  relative  error  in  the  same  stress  componets  computed  by  means  of  an  indirect 
(postprocessing)  technique  [13]  is  shown  in  Tab.  6.2. 

The  improvement  is  very  substantial.  The  postprocessing  technique  yields  stress 
values  which  are  not  sensitive  to  the  meshes  and  the  error  is  of  the  magnitud  lelgR. 


6.2.  Performance  of  the  p-Version 

When  the  solution  is  smooth,  the  p-version  performs  well.  When  the  solution  is 
not  smooth  and  the  elements  are  large,  then  the  pollution  error  is  generally  large. 
Satisfactory  theoretical  analysis  is  not  available.  It  is  known,  however,  that  the 
pollution  error  can  be  reduced  very  substantially  by  surrounding  points  of  stress 
singularity  with  one  or  more  layers  of  elements. 

The  relative  errors  at  point  Xj  »  0.1,  x2  “  0*2  (which  is  located  at  element 
boundaries)  are  shown  in  Tables  6.3a  and  6.3b.  The  results  presented  in  Table  6.3a 
are  strongly  affected  by  pollution  because  the  vertex  of  the  neighboring  element  was 
on  the  singular  point.  The  results  presented  in  Table  6.3b  are  much  less  affected  by 
pollution  because  an  extra  layer  of  elements  were  added  so  that  the  neighboring 
element  no  longer  had  a  vertex  on  the  singular  point.  The  local  error  is,  of  course, 
the  same  in  both  cases.  Table  6.3a  and  6.3b  illustrate  the  importance  of  proper  mesh 
design  when  the  stresses  are  computed  from  the  finite  element  solution  directly.  The 
postprocessing  method  removes  sensitivity  to  mesh  design. 
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TABLE  3.2a.-  H-EXTENSION  WITH  ASYMPTOTICALLY  OPTIMAL  MESH,  a  =  0.7 
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TABLE  3.2b.-  H-EXTENSION  WITH  ASYMPTOTICALLY  OPTIMAL  MESH,  a  -  1.1 
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TABLE  3.3.-  PERFORMANCE  OF  AN  OPTIMAL  MESH  (M  -  2,  a  -  1.1) 
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TABLE  3.4.-  MESH  PARAMETERS  (a  =  0.7) 
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TABLE  3.5.-  PENALTY  FOR  USING  OVERREFINED  MESH  IN  TERMS  OF  INCREASED  VALUES 
OF  C  (FIRST  ROW)  AND  N  (SECOND  ROW) 
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TABLE  3.6.-  P-EXTENSIONS  ON  UNIFORM  MESHES  (o  -  .7) 
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TABLE  3.7.-  PERFORMANCE  OF  p-EXTENSION  BASED  ON 

OPTIMAL  p— DISTRIBUTION  AND  UNIFORM  MESH 
(M  =  2) 
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3.8.-  PERFORMANCE  OF  h-p  EXTENSION  (a*  0.7,  s  =  0.4), 
ASYMPTOTICALLY  OPTIMAL  MESH  AND  p-DISTRIBUTION 
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TABLE  3.9.-  PERFORMANCE  OF  THE  h-p  EXTENSION,  ASYMPTOTICALLY 
OPTIMAL  MESH  AND  UNIFORM  p  (a-  .7) 
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TABLE  3.10a.-  H-EXTENSION  UTILIZING  FEEDBACK,  a  -  0.7 
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TABLE  3.10b.-  H-EXTENSION  UTILIZING  FEEDBACK,  a  *  1.1 
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TABLE  4.1.-  ERROR  IN  STRESSES  IN  DEPENDENCE  ON  POLYNOMIAL  DEGREE  AND  SMOOTHNESS 


8 


TABLE  5.4.-  P-EXTENSION  BASED  ON  PROPERLY 
DESIGNED  MESH 


p 

'“'er* 

.  1 

31.96 

2 

12.36 

3 

6.197 

4 

3.277 

5 

2.131 

6 

1.436 

7 

1.128 

8 

.8839  *  . 

TABLE  6.1.-  RELATIVE  ERROR  IN  STRESS  COMPONENTS  AT  POINT  (0.25,  0.25) 
DIRECT  COMPUTATION 


N 

■M 

„R  Z 
e22* 

-R  z 
e  12* 

A 

10.99 

7.57 

11.16 

1 

4.79 

2.32 

35.86 

221 

2 

12.21 

2.01 

.093 

3 

17.27 

17.42 

106.43 

4 

9.71 

13.01 

71.49 

A 

4.09 

5.47 

13.42 

1 

1.46 

2.68 

22.96 

617 

2 

4.41 

.099 

55.69 

3 

4.41 

13.43 

3.88 

4 

6.07 

11.74 

28.85 

TABLE  6.2.-  RELATIVE  ERROR  IN  THE  STRESS  COMPONENTS  AT  POINT 
(0.25,  0.25)  INDIRECT  COMPUTATION 
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TABLE  6.3a.- 


RELATIVE  ERROR  IN  STRESSES.  A  VERTEX  OF  THE  ADJACENT 
ELEMENT  LIES  ON  THE  CRACK  TIP 


TABLE  6.3b.-  RELATIVE  ERROR  IN  STRESSES.  THE  ADJACENT  ELEMENT  IS 

SEPARATED  FROM  THE  CRACK  TIP  BY  ONE  LAYER  OF  ELEMENTS 


p 

-R  » 
e  11* 

_R  y 
e22* 

eR  T 
e12* 

1 

13.57 

33.108 

17.39 

2 

2.124 

6.976 

.6688 

3 

1.091 

3.843 

3.620 

4 

.9997 

1.923 

1.852 

5 

.2653 

.7836 

.7166 

6 

.1702 

.3928 

.3503 

7 

.0784 

.2123 

.1529 

8 

.04005 

.1256 

.07431 

